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We discuss inhomogeneous ground states of the Nambu-Jona-Lasino (NJL) and quark-meson 
(QM) model within mean-field approximation and their possible existence in the respective 
phase diagrams. For this purpose we focus on lower-dimensional modulations and point out 
that known solutions in the 2-fl and l-fl dimensional (chiral) Gross-Neveu (GN) model 
can be lifted to the to the 3-1-1 dimensional NJL model. This is worked out in detail for 
one-dimensional modulations and numerical results for the phase diagrams are presented. 
Focus is put on the critical point and on vanishing temperatures. As an interesting result 
the first order transition line in the phase diagram of homogeneous phases gets replaced by 
an inhomogeneous phase which is bordered by two second order transition lines. 



I. INTRODUCTION 

Until today the phase diagram of quantum chromodynamics (QCD) is subject to intense theoretical 
and experimental investigations (for dedicated reviews see Ref. Experimentally its structure 
at finite temperatures is explored by heavy ion collisions, which currently focus on the formation 
and properties of a strongly interacting plasma at large temperatures as well as on the search for a 
chiral critical point in the phase diagram. Theoretically ab initio calculations are limited to small 
net-baryon densities and as a consequence possible scenarios at non-vansihing densities and strong 
coupling are often discussed within phenomenological models. 

Since NJL-type models^ share global symmetries and the phenomenon of chiral symmetry breaking 
with QCD, they are - typically in mean- field approximation - widely used to study the phase 
diagram at moderate temperatures and densities. In this context e.g. color-superconductivity [l], 0] 
as well as the location of the critical point (s) have been addressed [3|]. In recent years this branch 
of models has been extended to Polyakov-NJL-type models, which include the dynamics of the 
Polyakov-loop as the order parameter of confinement 0, 0]. However for vanishing temperatures 
these models reduce to the respective NJL-type model. 

In this work we investigate inhomogeneous ground states in the NJL Q] and QM model 0, Q]. 
These are characterized by a spatially varying order parameter and have been discussed for QCD 
at least in the large limit, where they are expected to form the ground state at sufficiently high 
densities 0,0, 
matter picture 



Related to this they show up in holographic models [12] and in the quarkyonic 
that suggests a similar structure for QCD. The investigations of these phases is 
however limited, mainly because they are technically much more involved. Within the NJL model 
such phases have been analyzed at vanishing temperatures applying further truncations [3] as well 



^ We refer to NJL-type models as models that at least in the applied approximation reduce to the NJL model, 
possibly extended by additional point-like interactions, on a technical level. This includes e.g. simplified ansatze 
for the gluon interaction, the use of different regularizations, the instanton liquid model and the QM model. 
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as for the so-called chiral density wave [l^j 0) • latter case the order parameter is assumed 

to be a plane wave and can be solved on mean- field level for vanishing current quark masses. In 
addition it has been suggested recently that the first order phase transition in the phase diagram 
of the NJL model, at least in the vicinity of the chiral critical point, is replaced by two second 
order phase transition lines that border an inhomogeneous phase (171] . The intersection of the two 
lines defines a Lifshitz point, which coincides with the critical point in the NJL model. 
Inhomogeneous phases have also been investigated in the 1+1 dimensional (chiral) GN modeP 
for the large limit, which technically corresponds to a mean-field approximation. Not only the 
understanding of the phase diagram has been improved significantly as it possesses regions with 
inhomogeneous ground-states 18|, Il9l . |20|], also theoretical understanding of the integrability of 
these models has been obtained. In the chiral limit all self-consistent solutions can in principle be 
classified [2l| and the phase diagrams have been analyzed in detail 22]. For the case of the GN 
model it is furthermore possible to introduce finite quark masses and to study their effect on the 
structure of the phase diagram [3] • 

The underlying idea of the present investigation is based on the observation that the theoretical 
problem of finding self-consistent inhomogeneous phases with a lower dimensional modulation can 
be reduced to a problem in a lower dimensional model. This is mainly due to the structure 
of the mean-field Hamiltonian, which is of Dirac-type. As a consequence self-consistent phases 
with a one-dimensional modulation in the NJL model can be studied on similar grounds as in 
the 1+1 dimensional (chiral) GN model and all results obtained for the latter can be used here. 
In particular all self-consistent inhomogeneous solutions with a one-dimensional modulation are in 
principle also known for the NJL model, at least in the chiral limit. Furthermore we can investigate 
self-consistent inhomogeneous phases at finite quark masses, which has not been possible before. 
The considered inhomogeneous ground states are then lattices of domain- wall solitons. In the same 
way self-consistent phases with a two-dimensional modulation could be related to those in the 2+1 
dimensional GN model. However no analytical solutions for a real two-dimensional modulation are 
known for the latter. 

Being able to investigate the role of inhomo gen eous phases in the phase diagram of the NJL model, 
we can confirm the picture obtained in Ref. [17| for the vicinity of the critical point and the absence 
of a first order phase transition line in the phase diagram^. Furthermore we can analyze the relation 
between inhomogeneous phases and e.g. the strength of the first order phase transition (present in 
the case of homogeneous phases) . For the main part of the paper we do not allow for pseudo-scalar 
condensates, especially since only without those self-consistent solutions for finite current quark 
masses are known. We will however also discuss why this condensates are not expected for the 
ground-state. 

At the end of the paper we will finally also consider the QM model. The purpose of this is twofold: 
On the one hand we would like to extend the analysis to a larger class of models in general, on 



^ The chiral GN model is also called the 1+1 dimensional NJL model. The global Z2 symmetry of the GN model is 

extended to a U{1) symmetry in this model. 
^ This picture could of course be modified by the inclusion of color-superconducting phases. 
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the other hand the regularization of the NJL model for inhomogeneous phases is non-trivial with 
regard to a combined vacuum and QCD phase diagram phenomenology. 

The paper is organized as follows: In section [TI] we introduce the NJL model and the applied mean- 
field approximation allowing for inhomogeneous phases. After this we show in section IIIII how 
self-consistent solutions for lower-dimensional modulations can be obtained from lower dimensional 
models. We then summarize self-consistent solutions of the (chiral) GN model in section IIVI and 
work out in section |V] how those can be used for the NJL model. Before coming to the QM model 
in section IVIIl we remark in section IVII why pseudo-scalar condensates are not expected. Finally 
we discuss various numerical results in section IVIIII and close with a summary and outlook in 
section IIXI 

II. NJL MODEL AND MEAN-FIELD APPROXIMATION 

In this work we concentrate on the two-flavor NJL model given by the Lagrangian 

£ = ^ {ij'^d^ - m) + ( (V^V^) ' + (^iT'r V) ') , (1) 

where ip is the 4A^jA'^c-dimensional quark spinor for Nf = 2 flavors and Nc = 3 colors, 'j^ are Dirac 
matrices, Gs is the scalar coupling and rh the mass matrix for degenerate quarks with current 
quark mass m. For A'^^ = 2 the matrices r° are the conventional Pauli matrices. 
We address phases with non- vanishing expectation values {ipip) = 5'(x) and {'tpi'y^T^ip) = Pa(x). 
Expanding bilinears around those expectation values and neglecting quadratic contributions, we 
then work within mean-field approximation. For technical reasons that will become more apparent 
in the following we furthermore restrict ourselves to the case where the direction of the vector 
Pa(x) is constant in flavor space so that we can choose a flavor basis where -Pi(x) = P2{x) = and 
P3(x) = P(x). The mean- field Lagrangian therefore takes the form 

Cmf = i^{il^d^ - m + 2G,5(x) + 2iG,7V='P(x))V' - (5(x)2 + P{^f) . (2) 

Since the mean-field Lagrangian is bilinear, the thermodynamic potential as an effective action 
in the expectation values can be formally evaluated. In the case of a periodic condensate with 
Wigner-Seitz cell V and using the imaginary-time formalism (see e.g. Refs. the thermodynamic 
potential as an effective action in ^(x), -P(x) is then given by 

0(r, S(x), P(x)) = In / Vi^V^lj exp f / {Cmf + ^^tV) | 

= E Tr^jy Log (i + Hmf " ^)) + ^ ^ ^ ^^"^'^ ' 

(3) 

where the functional logarithm and trace act on Dirac, color and coordinate space. The 
Hamiltonian Hmf is obtained from Eq.(l2]) and reads 

Hmf = -il^l'di + 7° (m - 2G,5(x) - 2iG,7V3p(x)) , (4) 
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which is a direct product of the two isospectral Hamiltonians 

Hmf,± = -ij^Ydi + 7° (m - 2G,5(x) ^ 2iG,75p(x)) . (5) 

Assuming a sensible regularization for the functional trace we can then express the thermodynamic 
potential through the eigenvalues {En} of Hmf,+ as 

n{T, ^; M(x)) = ^ Yl 1" (2 cosh iy^^j J^vJv 4^ ^ ' 

where M(x) = m — 2Gs{S{x.) + iP(x)). We note that the mass function M(x) as the order 
parameter can be complex and parameterizes a U{1) subgroup. 

In order to minimize the thermodynamic potential in the order parameter M(x) we can address 
the stationary constraint ^j^^^^t = 0, which in turn gives the gap-equation 

M(x) = „ + i^^ta„h(^)*t(x)|^A,(x) 

En 

as a self-consistency condition on M(x). Here ipni^) are the normalized eigenvectors of the Hamil- 
tonian for the eigenvalues En, i.e. y fyiptipn = 1. 



III. LOWER-DIMENSIONAL MODULATIONS 

Having the formal expressions for thermodynamic potential and gap-equation, we turn towards 
inhomogeneous phases. The focus is the determination of the eigensystem of a Hamiltonian and the 
check of self-consistency. As will be discussed in the following, these problems can be dimensionally 
reduced for lower-dimensional modulations by help of Lorentz symmetry. 

Suppose we have a representation S'(A) for the Lorentz transformation A that acts on the momen- 
tum operator of free spinors = {H, Pi)^ as 

S{A)P^'S-\A) = A'^^P'^. (8) 

If the system is translationally invariant in one or more directions, the corresponding momenta 
Pi commute with the Hamiltonian H and we can label its eigenstates also by p_L, the eigenvalue 
of Pj_. Now take 'i/'A,o to be the eigenvector with Hipx^o = Xip\^o and Pj_^a,o = 0, A^,y to be the 
Lorentz transformation that boosts (A, O)'^ to (^X\J 1 + Pj_ f^'^, and define 

V'AVTTiiTA^.o, = fVl + Pi/A^l '5-^(A)^A,o. (9) 



Px 



The prefactor implements the proper normalization for the conventional choice of 5(A) 23] and it 
is straightforward to check that 
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Therefore the whole eigenvalue spectrum can be constructed from the subspace spanned by {V'A.o} 
and for the thermodynamic potential we obtain 



17(r,^;M(x)) 



2TNc 
Vn 



E 



(27r)°'j 



1 f |M(x) - m\ 



In I 2 cosh 



+ const. , 



AJl + p2/A2-/. 



2T 



(11) 



where V|j is a Wigner-Seitz cell of the lower dimensional modulation. Addressing the self-consistency 
condition in Eq.© we note that 5(x) and -P(x) are scalar and pseudo-scalar, respectively. We find 

A 



■ ' \ • P 1' ^AVl+pi/A2 
^AVl+pl/A^Px^^'^AVl+pi/A^P 

and the gap-equation takes the form 



A 

A^pi 



= V'A,07 V'A.O 



(12) 



M(x) = m + 



(27r) 



tanh 



AJl + pi/A2-^ 



2T 



x^ + pI 



-4\,o (1 - 7^) V'A.O • 



(13) 



Therefore we have formal expressions for the thermodynamic potential and the gap-equation in 
terms of the subspace spanned by the eigensystem corresponding to p_L = 0. To be more specific 
we now choose the Weyl representation for the 7-matrices, i.e. 
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1 

1 



7 



a' 
-a' 
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-1 
1 



(14) 



with the conventional Pauli matrices cr* so that the Hamiltonian in Eq.([5]) takes the form 



Hmf,- 



ia'di M(x) 
M(x)* -ia^di 



(15) 



First restricting to a one-dimensional modulation in the z-direction, the Hamiltonian in the sub- 
space corresponding to p_L = is 

/ id, M{z) \ 

-id, M{z) 
M{z)* -idz 

M{z)* id, j 



HmF: 



ID 



which can be block-diagonalized into 



Hi 



MF-AD 



\ 



^ Hid(M{z)) 



(16) 



HiD{M{zy) 



(17) 
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(18) 



where 

\ M{z)* id^ J 

The problem is therefore reduced to finding the quasi-particle spectrum of Hid which is the Hamil- 
ton of the chiral GN model and which will be discussed in more detail in section [TVl 
Considering the scenario with a two dimensional modulation it is more convenient to assume the 
system to be translationally invariant in the z-direction. The Hamiltonian in Eq. ()15p for = 
is then given by 

/ id,-dy M{x,y) \ 



HmF:2D 



idx + dy 
M{x,y)* 



\ 



M{x,y)* -id^-dy 



M{x,y) 
-idx + dy 



(19) 



This can be block-diagonalized only for real M(x,y), i.e. for P(x) 
the problem into the form 



0. In this case we can cast 



/ -M{x,y) idx + dy 



HmF:2D 



\ 



idx-dy -M{x,y) 



\ 



(20) 



M{x,y) idx + dy 
idx - dy M{x, y) j 

which corresponds to the GN model in two-dimensions. The fact that we need a real mass function 
M{x, y) is not surprising since there is no chiral GN model with two-component spinors in two 
dimensions. 

In general we are therefore able to constrain ourselves to a lower dimensional model if we want 
to determine the quasi-particle spectrum for an inhomogeneous phase with a lower dimensional 
modulation. However a self-consistent modulation in the lower dimensional model does not neces- 
sarily have to be self-consistent in the higher dimensional model. For the cases considered in this 
paper and essentially all analytical solutions known this will nevertheless be the case since for each 
eigenvector separately we will have 



V'A.o (l - 7^) V'A.o oc a + bM{z) 



(21) 



and there are enough parameters in the analytical expression for the modulation to tune the 
coefficients a, b such that the solution is self-consistent. 

IV. INHOMOGENEOUS PHASES IN ONE-DIMENSIONAL MODELS 



Having realized that self-consistent solutions with lower dimensional modulations may be obtained 
from lower dimensional models, we now present some self-consistent phases with a one-dimensional 
modulation. This reduces to the study of the (chiral) GN model whose Hamiltonian is given by 
Eq. pSj) and in the following we summarize some results deduced from Refs. [isl. [l9. 
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The thermodynamic potential can be evaluated by the knowledge of the density of states p(A) 
which enables us to perform the sum over the eigenvalue spectrum: For an arbitrary function /(A) 
it is given by 

f = [dXpWfiX). (22) 

II X 

Recent investigations have shown that at least for the chiral limit all self-consistent solutions can 
be classified and that all information about the inhomogeneous phase is encoded in the band-edge 
energies {Ei} of the valence bands of the quasi-particles For the case of one gapped valence 
band the general solution for a complex order parameter M{z) has been determined 2M the 
spectral density associated with the Hamiltonian Hi£) is given by 



J_ a{\)+{\M{z)\^) jj4 
21- /Tfi i\_T^:\ ' j 



m^^{\-Ej)>Q 



p(A) = { ^^U^-^^) ' ' , (23) 

, else 



where 



4 

a(A) = 2(A-iX:^.y-^E(^^-^^)' 



j=l i<j 

1 



(El - E^)iE2 - EAm %^~f\\f-f\ ) 

'^y{Ei-Es)(E2~E4)) 

and El < E2 < E^ < Ea are the edges of the valence bands^. For more details on the function 



M(z) we refer to Refs. 2^, l2J] since we will limit ourselves to real order parameters M(z) in the 
following. 

For real order parameters M(z) we have p(A) = p(— A) or equivalently Ei = —E4, E2 = —E3. 
So there are only two parameters to specify the solution, which we choose to be given through 
£"3 = \/l — i/A and E4 = A [3]. This case is of particular interest since it is possible to generalize 
the solution to finite current quark masses [3]. The order parameter then has the specific form 

M{z) = Afz.sn(6|z.)sn(Az|zy)sn(Az + 6|zy) + ^^^^^4j?^^^P^^ > (25) 
V sn(6|i/) ; 



where A is a scale parameter and sn, cn, dn are elliptic Jacobi functions with elliptic modulus y^, 
which physically describes lattices of equidistant solitions. Due to periodicity the parameters b and 
v can be limited to 6 G [0,K(i/)] and G [0, 1] with the quarter period K(z^). The chiral limit, i.e. 
vanishing quark masses, corresponds to b = K(z/). 

In terms of the eigenvalue spectrum, the case of finite current quark masses included through 
b 7^ K(i^) only leads to a shift in the distribution of the eigenvalues A 



A ^ sign(A)\/A2 + <5A2, 

S = 2m ^ - 1 ^ (26) 



The expression for (|Af(z)p) has not been derived in Ref. [2J], but is obtained by solving the Ginburg-Landau 
equations and presenting the solution in terms of Jacobi elliptic functions. Since it is not relevant for the main 
focus of this work, we do not give the tedious derivation. 
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and as a consequence the sum over the eigenvahie spectrum for arbitrary functions /(A) can be 
performed by using the density of states p(A) of the chiral hmit 



A 



1 /""^ 

= / d\p{\)f{X) 

II \ J —oo 

/oo 
dX p(A)/(AVl + 5A2/A2) . (27) 
-oo 

Addressing self-consistency and following Ref. [3], we have for the normalized eigenvectors i/jioiz) 
of Hid 

J . . , A AV^((^ + l)(J-i. + l) 

Mz)^Mz) = - SA^ - A2E(.)/K(.) '''^'^ - X{X^-6A^-A^B{u)/K{u)) ' ^'^^ 

As discussed in the context of Eq. (j2ip . we observe that the first term in Eq. (j28p is proportional 
to the modulation M{z) whereas the second is constant for each eigenvector. This enables us to 
obtain self-consistency in the gap-equation for a proper choice of parameters. In particular we find 
that the second term vanishes for 6 = which then corresponds to the chiral limit. Furthermore 
we have 

V^A.o (1 - 7^) ^A,o = V^idV'id (29) 

and conclude with the same reasoning that we can achieve self-consistency for M{z) being a one- 
dimensional modulation in the NJL model. 



V. ONE-DIMENSIONAL MODULATIONS IN THE NJL MODEL 

We will now concentrate on real order parameters M(z) and parametrize the band edges by = 
\/l — uA, = A, which corresponds to the order parameter in Eq. (j25p for b = K(z^). After a 
non-trivial computation using Eq. ()23p we obtain the density of states 

r 1 y-A^EH/K(.) , (A2 - A2) (A2 - (1 - z.)A2) > 

I , else 

which of course agrees with the finding in Ref. The evaluation of the thermodynamic potential 
in Eq. (jlip then reduces to evaluating expressions of the form 



^ clEE^ j dup{Eu)f (^^/WTdA^^ 

dE ~p{E)f (yW^^i^^ . (31) 
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Here we introduced f{x) = f{x) + /(— x), used p{X) = p(— A) and implicitly defined the effective 
density of states 

1 

p{E) = — / duE^p{Eu). (32) 

^TT J_i 

The overall factor of two in Eq. (j3ip stems from the degeneracy of the two Hamiltonians Hio{M) 
and Hi£){M*) in Hmp-h). A straightforward but somewhat tedious calculation using the density 
of states in Eg . ([30]) then gives the explicit expression 



_ B{e\9) + EH/KHF(g|z>) - F{9\9) + y^{E'^ - A'^){E'^ - 9A'^)/{EA) ^ 

(33) 

where we introduced the elliptic integrals of 1st and 2nd kind as well asz> = 1 — v, 9 = 
arcsin(i?/(\/^A)) and 6 = arcsin(A/£'). As an interesting cross-check we can consider the limit 

mi=i= ^9{E-A)^{E^-A^)E, (34) 
which is the effective density of states in an homogeneous phase with a quasiparticle gap A and 

mUo = (35) 

which corresponds to the ultra-relativistic gas. We also find for the asymptotic behavior of the 
effective density of states 

P^^^ - -^-^^ + )■ 



For the considered one-dimensional modulation the thermodynamic potential in Eq. ijlip can there- 
fore be cast into 



\z] — mP 



{T,p;A,u,S) = -2N, dE p{E) f^^,, (^/ E^ + 5A^) + ^ dz\M{: 

+const. , (37) 



where L is the period of the modulation and 



/bare(2;) = T hi ( 2 cosh ( ^ ) ) + T In ( 2 cosh ^ 



2T J J \ \ 2T 

= fvv{x) + /medium (a;) , 

/uv(2;) = X, 

Lcdiumix) = Tln(l + exp(-^))+Tln(l + exp(-^] ] . (38) 
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The last missing step is then a regularization of the diverging integration. For homogeneous phases 
this is mostly done by a momentum regularization (see e.g. Refs. Q, 0]). This is however not 
possible for inhomogeneous phases since the quasi-particle energies can no longer be labelled by a 
conserved three-momentum. Instead we have to apply a regularization of the functional logarithm, 
e.g. by a proper-time regularization, which is essentially a regularization acting on the energy 
spectrum instead of the quasi-particle momenta. For this purpose we have already identified the 
divergent vacuum contribution associated to f\jY{x), which we regularize by a specific blocking 
function in the proper-time integral leading to a Pauli-Villars regularization of the form [25 1 

3 

fvvix) ^ /pv(x) = ^CjVx^+j^\ (39) 

i=o 

with Co = 1, ci = —3, C2 = 3, C3 = —1 and a cutoff scale A. 

With the expression of the thermodynamic potential as a function of A, and 6 we are now able to 
perform numerical investigations by performing the energy integration numerically and minimizing 
in the mentioned parameters. The resulting phase diagrams will of course depend on the choice of 
the model parameters Gs, A and m. Those can be related to the chiral condensate {tpip) and the 
pion decay constant though the expressions 25 1 

j=0 \ 1 J 

where M„ is the constituent mass for the chirally broken phase in the vacuum. 




VI. A NOTE ON PSEUDO-SCALAR CONDENSATES 



Having worked out an expression for the thermodynamic potential for real order parameters -/Vf(z), 
we may also discuss the more general case of a complex order parameter or in other words the 
relevance of pseudo-scalar condensates Paix). 

On the one hand side this is less attractive from a technical point of view, since we only know self- 
consistent solutions away from the chiral limit for real order parameters M(z). More importantly 
we can also argue that these phases are energetically less preferred, at least in the vicinity of a 
second order phase transition to the chirally restored phase and in the chiral limit. 
For this purpose we consider the generalized Ginzburg-Landau (GL) functional which is a system- 
atic expansion of the thermodynamic potential in the magnitude of the order parameter M(z) as 
well as in gradients acting on it. Both are treated to be of the same order. For the chiral GN 



11 



model this has been worked out in Ref. 26|], where 



nGL,GN{M) - nGL,GN{0) = -j\Mf + -^lm{MM'*) + -j{\M\^ + \M' 

+^Im((M" - 3|MpM)M'*) 

+—(\Mf + SImPIM'P + 2|mP|m2|' + -|M"P) + . . . (41) 
6 2 

has been derived. The dependence on temperature and chemical potential is hidden in the GL 
coefficients Oj and the details are irrelevant for the present discussion. For the NJL model the GL 



expansion for real order parameters has been worked out in Ref. [17| and a remarkable similarity 
to the GN model has been found when restricting to one-dimensional modulations. From the 
viewpoint of Lorentz symmetry as presented in section Hill this may not be surprising since the GL 
coefficients are defined as integrals over functions in momentum space, so that at least for = 
the integrands of chiral GN and NJL model have to be equal up to an overall prefactor. Since the 
rotationally invariant extension to p_L 7^ is unique, the structure of the GL expansions has to be 
the same. 

In contrast to the chiral GN model the NJL Hamiltonian is however the direct product of Hio[M) 
and Hi£,[M*). Therefore the structure of the GL functional for the NJL model limited to one- 
dimensional modulations corresponds to the sum of that in Eq. (j4ip and its conjugate. We therefore 
obtain 

nGL,NJL{M) - nGL,NJL{<^) = y + ^(|M|^ + |Mf) 

+—{\Mf + 3|Mp|M'p + 2|Mp|MY + -\M"\^) + . . . , (42) 
6 2 



where (3i are the GL coefficients of the NJL model as stated in Ref. 17l |. 

With the generic form of the GL functional at hand, we can now compare various ansatze for M{z). 
We will limit ourselves to second order phase transitions from an inhomogeneous to the chirally 
restored phase and therefore to a regime where /32 > 0, /34 < 0. Furthermore all Fourier modes 
of M(z) decouple to order due to momentum conservation in the GL coefficients, which are 
calculated in an homogeneous background. For this reason we consider only the two extremes 

Mff{z) = ^^ey.v{iqz), 

Msin{z) = Asj„sin(g2;). (43) 

For an individually optimized wave- vector q the thermodynamic potentials for these ansatze is then 
given by 

^gl,njl{Mff)-Qgl,njl{^) = A2,^-i/34A^^ + 0(Af.^), 

^GL,NJL{Ms^n) " f^GL,JVJL(0) = (f " ^) ^In " Y^I^AK.n + 0(AL) • (44) 

From this result we can deduce that the phase transition for both cases happens at the same 
values of the GL coefficients, namely at /32 — = 0. The inhomogeneous phase then occurs for 
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h — 1^ < 0, where the coefficient of the ^\p- and A^.^-term is negative. Since the coefficient 
of the A^j„-term is however smaUer than that of the A^^-term (recah that /54 < 0), we find that 
the modulation of the form Msin{z) is energetically preferred compared to Mff{z) for optimized 
parameters. In this sense the ground state with a real order parameters wins over the plane wave. 
Two brief comments are in order here: First, the self-consistent modulations discussed in the 
previous section indeed become sinusoidal when approaching the phase transition to the chirally 
restored phase and second, the sinusoidal modulation can of course be understood as two plane 
waves in opposite directions. The latter have also been found energetically preferred compared to 
the plane wave in inhomogeneous (color-)superconductors [271. l28l|. 



VII. QM MODEL AND MEAN-FIELD APPROXIMATION 

The NJL model regularized by a proper-time regularization and adjusted to chiral condensate and 
pion decay constant is known to give constituent quark masses of order 200MeV in the vacuum. 
Hence it gives an undesired phenomenology with regard to the QCD phase diagram, mainly be- 
cause quasi-particles will start forming a Fermi surface at ^ ~ Mq. This is phenomenologically 
unacceptable for ^ < (Mtv — i?)/3 ~ 308MeV, where M^v is the nucleon mass and B the binding 
energy of nucleons in nuclear matter. In contrast the NJL model regularized by a sharp three- 
momentum cutoff gives constituent quark masses in the vacuum significantly larger than 300MeV. 
Since we haven't found a regularization that avoids these problems, we instead also introduce a 
model that is very similar to the NJL model and where the issue of the regularization scheme can 
be surpassed: The linear sigma model, which in this context is usually named QM model 29]. 
The Lagrangian of the QM model with Nf = 2 and Nc = 2>\s given by 

Cqm = ^ {il^d^ - g{a + ^75^^^)) - U{a, vr'^) , (45) 

where 

U{a,Ti'') = -i(a^aaV + a^7r"a^^") + ^(a2 + 7rV-z;2)'-ca, (46) 

tjj is again the 4A'^jA''c-dimensional quark spinor, a the scalar field of the u-meson and vr" the 
pseudo-scalar fields of the pion triplet. In mean- field approximation we treat the fields a and vr*^ as 
classical and replace them by there expectation values [a, l29j. Furthermore we can use low-energy 
relations to connect the parameters c, g, X and with hadronic observables. We will express those 
by the pion-decay constant /tt, the constituent quark mass in the vacuum Mq, the pion mass m.^ 
and (7-meson mass ttIo- via {a) = {it"') = 0, c = m^fT^, g = Mq/ f^^, A = (m^ — m^)^/(2/^) and 
v"^ = fl- ml/\. 

For the thermodynamic potential in mean-field approximation we only include the contributions 



13 



of the fermionic fluctuations and approximate 

0QAf(T,^;a(x),7r"(x)) = -^In [ V^VipVaVn^ exp I [ (£QAf + /^V^tV) 

^ J \Jxe\o,l]xV 



-field TNc 



V 

n 



TtDJ.V Log i ^ [iiOn + Hmf,qm - 



+ i |^C/(a(x),vr'^(x)), (47) 

where cr(x) and 7r'^(x) are non- vanishing expectation values of the respective fields. Limiting to 
cases with '7t^{'x.) = vr^(x) = the Hamiltonian reads 

Hmf,qm = -ilSdi + 7° (5^(x) + i57'T=^7r3(x)) , (48) 

and we have to evaluate the same functional trace-logarithm as in the case of the NJL model, but 
now with the identification Af (x) = gf((T(x) + i7r^(x)). 

The QM model is renormalizable which means that the divergences in the functional trace- 
logarithm can be absorbed by the model parameters. For this purpose it is useful to observe that 
we can always separate zero-point fluctuations from thermal fluctuations as indicated by f\]\{x) 



and /medium C 


x) in 


follow Refs. [ 


^1 


29 



29|], where it has been assumed that the contribution from zero-point approxima- 
tions can well be approximated by y Jy [/(ct(x), 7r'^(x)) with the parameters directly adopted to 
phenomenology. As a result we can evaluate the thermodynamic potential for one-dimensional 
modulations cr(x) = M{z)/g and 7r"(x) = with M(z) given in Eq. ()25p on the same level as for 
the NJL model and we get 

^MF,QMiT,n;A,u,6) = ^}MF,QM{T,fj,;M{z)/g,0) 

= -2N,j'" dEpiE) U,dinm (Ve^ + 5A^) +jj'^dzU{Miz)/g,0) 
+const. . (49) 

For completeness we have to denote that this will only be a self-consistent solution if M(z) would 
solve the equations of motions associated to minimizing the functional [/(^(x), 7r"(x)), which is 
in general not the case. This is not considered a major issue within this work since we already 
limited the space of possible solutions to that of one-dimensional modulations so that we cannot 
guarantee to find the true ground state anyway. Instead we will find domains where at least one 
inhomogeneous phase is energetically preferred over all homogeneous phases and therefore pointing 
out domains where the ground-state is for sure inhomogeneous (in the presented approximation). 
At the end of this section we like to point out that we can easily analyze the vicinity of the chiral 
critical point by a generalized Ginzburg-Landau approximation. For this purpose we first observe 
that the functional trace-logarithm in Eq. (j47p is formally identical to the the one in Eq.(l3|). The 
GL expansion of the latter has been worked out in Ref. [171] and the structure of the GL functional 
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stays the same when limiting to thermal fluctuations^. The GL expansion of ^mf,qm{T, /i; <7(x) = 
M{x.)/g, 0) can therefore easily deduced from Ref. [l7[. Without giving unnecessary details for the 
presented discussion, the GL expansion with coefficients (5[ takes the form 

nMF,QM{T, a(x) = M(x)/5, 0) = ^J^[ f M(x)2 + ^(M(x)4 + (VAf(x))2) 

+^(M(x)6 + 5M(x)2(VM(x))2 + (AM(x))2) 



1 



cgM{x) - ^L^M{x 



,2 



+^{g'XMi^f + 252(VM(x))2)j + . . . . (50) 

Restricting to the chiral limit for simplicity, i.e. c = 0, we can identify the chiral critical point as 
the position where the coefhcients of the M(x)^-term and the M(x)^-term vanish, i.e. 

/3^-5VA = 0, 
Pi + g^X = 0. (51) 

Considering the coefficient of the (VM(x))-^-term at this point we find 

being negative for A > 1. Since homogeneous phases are instable for a vanishing M(x)^-term 
and a negative (VM(x))^-term, the chiral critical point would be hidden in this case since an 
inhomogeneous phase were energetically preferred at its location. This scenario is actually real- 
ized for phenomenological parameters since we have A = m^/(2/^) ~ 20. Due to the involved 
approximations we do however not want to overemphasize this result. 



VIII. NUMERICAL RESULTS 

With the expressions for the thermodynamic potentials of N JL and QM model in Eq. (I37p (combined 
with the described regularization) and Eq. (j49p at hand, we can now study the structure of the phase 
diagram numerically. Except for the analytically known averages of M{z), M{z)'^ and M'(z)^ given 
in Ref. [3], we evaluate all involved integrals numerically and obtain nj^.[pj^jiiQj^,.j{T^^;/S.,v^5). 
Since numerical minimizations typically determine parameters up to a total and/or relative error 
and since elliptic integrals depend on the elliptic modulus in a highly nonlinear fashion for 
— > 1, it is numerically very helpful to express v = 1 — 16exp(— 41n(2)/q/i) with qh £ [0, 1] and 
then minimize in the parameters A, qf^ and 6. 

We will now start discussing results for the NJL phase diagram in the chiral limit, i.e. for m = 0. 
For the reasons discussed at the beginning of section IVIH we will fix = 88MeV to its value in 



^ The GL coefficients are defined through integrals in momentum space for which we can separate vacuum and 
thermal contributions. 
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FIG. 1: Left: Structure of the NJL phase diagram in the chiral hmit as a function of temperature T and 
quark chemical potential fiq for Mq — SOOMeV. The black (short-dashed) line indicates the second order 
phase transition from chirally broken to restored phase, the red (solid) line the first order phase transition 
and the bullet the critical point. The spinodal region is enclosed by the blue (long-dashed) lines. Right: 
Same plot as on the left including the orange (shaded) domain where the energetically preferred ground 
state is inhomogeneous. 



the chiral limit and instead of adjusting {ipip), we choose a value for Mq. These choices in turn fix 
the model parameters through Eq. (j40p . On the left side of Fig. [T]we present the phase diagram 
restricting to homogeneous phases for a value of Mg = SOOMeV. For larger temperatures we see the 
second order phase transition line from the chirally broken to the restored phase, which turns into a 
first order phase transition line at (^Ucr = 269MeV, Tcr = 74MeV) locating the chiral critical point. 
The first order line then ends at (/j, = 312MeV,T = OMeV), where the associated spinodal region 
spannes about 18MeV in the quark chemical potential fiq. Although the chiral condensate with 
{^/j^) = -(193MeV)3 for this case is phenomenologically too small, we observe that the structure of 
the phase diagram is similar to that typically found in NJL models using a sharp three-momentum 
cutoff 3. 

Focusing on the region near/around the first order phase transition and the critical point, the right 
hand side of Fig. [T] shows the same lines as on the left, but now also including the domain where 
inhomogeneous phases are energetically preferred. As discussed in Ref. 17|] for the vicinity of the 
critical point, we observe that there is no longer a first order phase transition in the phase diagram, 
since it is replaced by an inhomogeneous ground state. The transitions from the chirally broken to 
the inhomogeneous and from the inhomogeneous to the restored phase are both second order, where 
the first transition is characterized by the formation of (in the perpendicular direction) localized 
domain-wall solitons and the second by the melting of the condensate. 

The nature of the phase transitions is also apparent in the squared spatial average (M(a;)^)2 
of the order parameter as well as the wave- vector of the one-dimensional modulation q. For 
vanishing temperatures these are depicted on the left of Fig. [2] for Mq = 250, 300, 350MeV. We 
observe that q continuously raises from q = at the transition from the chirally broken to the 
inhomogeneous phase, which is related to the formation of localized objects. On the other hand 
{M{x)'^)^ continuously goes to zero at the transition from the inhomogeneous to the chirally 
restored phase. In the same plot we also see that at vanishing temperatures the constituent quark 
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FIG. 2: Left: Wave-vector q (dashed lines) and average of constituent mass -\/ {M{z)^) (solid lines) at 
vanishing temperatures as function of quark chemical potential /i^ for Mq — 250MeV (black lines), Mq = 
SOOMeV (red lines) and Alq = 350MeV (orange lines). Right: Same plot as on the right of Fig. [U now 
including results for Mq = 350MeV (upper branch) and Mq = 250MeV (lower branch). 



mass stays constant for < Mq (below any phase transition) and starts decreasing with the 
condensation of quarks. 

Whether or not quarks condense before reaching a phase transition has tremendous consequences 
on the structure of the phase diagram as shown on the right hand side of Fig. [2l Here we present 
the region of the first order phase transition for Mq = 250, 300, 350MeV. Qualitatively the picture 
for those three cases is similar, but quantitatively there are significant differences. Since for Mq = 
250MeV the quarks already form a Fermi surface before reaching the phase transition, the order 
parameter is already diminished at the phase transition. As a result the phase transition is very 
weak and the spinodal region/domain of inhomogeneous phases strongly reduced. For Mq = 
350MeV on the other side the interaction is very strong and there is no condensation before reaching 
the onset of inhomogeneous phase (or the first order phase transition when limiting to homogeneous 
phases). The domain of inhomogeneous phases is very large here and it is conceivable that its 
persistence beyond ^q = 400MeV is a regularization artefact. We also point out that although the 
value of constituent quark mass Mq is comparable to that in NJL models using a three-momentum 
cutoff, the chiral condensate with (tp^) = — (186MeV)^ at Mq = 350MeV decreases with increasing 
Mq here. 

The results indicate that the size of the domain for inhomogeneous ground-states, similar to the 
spinodal region, is related to the strength of the first order phase transition. In order to give a better 
picture for the latter, the thermodynamic potentials as effective actions in the order parameter 
are shown on the left of Fig. [3] at vanishing temperatures and the respective first order phase 
transition points. We find that the local maxima between Mq = 250MeV and Mq = 350MeV in the 
potential even when measured in the respective bag constants B^ac = ^(0, 0; 0) — 0(0, 0; Mq) differs 
by more than an order of magnitude, therefore underlining the quantitative difference between 
the two cases. For completeness, the value of the bag constant for Mq = 250, 300, 350MeV is 
B^ac = 40,62,91MeV/fm3, respectively 

Coming back to Mq = 300MeV we also like to give an impression of the thermodynamic poten- 
tial at vanishing temperatures using the right hand side of Fig. [3l Here we show the value of 
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FIG. 3: Left: Difference of thermodynamic potential Vl{T,ii]M) in the chiral hmit for the homogeneous 
phase with constituent mass M and its value for M = in units of the bag constant B^ac- The black 
(long-dashcd) line corresponds to Mq — 250MeV, the red (short-dashed) to Mq — SOOMeV and the orange 
(solid) to Mq — 350MeV. All for vanishing temperatures and /ig at the first order phase transition (when 
limiting to homogeneous phases). Right: Difference of thermodynamic potential ^phase in the chiral limit 
for various phases and its value for M = for Mq = SOOMeV in units of the bag constant Byac- The black 
(long-dashed) line corresponds to the inhomogeneous, the red (short-dashed) to the chirally broken and the 
orange (solid) to the chirally restored phase. 
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FIG. 4: Left: Same plot as on the right of Fig. [l] now including the domain of inhomogeneous phases 
for m = 5MeV and m ~ lOMeV. Branches with critical points at smaller temperature T and larger quark 
chemical potential /iq correspond to larger current quark masses m. Right: Same plot as on the left for 
Mq = 350MeV. 



the thermodynamic potential for the energetically most preferred homogeneous phase as well as 
inhomogeneous phase. Plotted are only cases that form a local minimum of the thermodynamic 
potential as an effective action, such that e.g. the spinodal region can be deduced from the plot. 
In the case of several local minima at a given the global minimum is energetically preferred and 
we see that if existent, the inhomogeneous phase forms the ground state. Also the order of the 
phase transitions can be deduced from the plot. 

As a next point we want to address the case of finite current quark masses. For this purpose we 
again fix the parameters Gs-, A by choosing /jr, Mq in the chiral limit and then turn on a finite 
current quark mass m. In Fig. |3]we present our results for the relevant part of the phase diagram 
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FIG. 5: Left: Phase diagram for the QM model in the chiral hmit with the red (dashed) hne indicating the 
first order phase transition when hmiting to homogeneous phases only. Inhomogeneous phases are preferred 
in the orange (shaded) domain enclosed by black (solid) lines. Right: Same plot as on the left, now for 
m„ = 69MeV. 
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FIG. 6: Same plot as on the left of Fig.[5l now for = 138MeV and with the bullet showing the critical 
point. 



choosing m = 0, 5, lOMeV. Aside from the fact that the second order phase transition is turned 
into a cross-over when leaving the chiral limit, we see that the qualitative picture stays the same: 
Again the first order phase transition line is replaced by two second order phase transition lines 
which enclose a domain where inhomogeneous phases are energetically preferred. All lines meet at 
the critical point, which can be calculated on the bases of homogeneous phases only and which is 
known to shift towards smaller temperatures and larger quark chemical potentials when increasing 
m. 

At the end of this section we sketch some results for the QM model and due to the similarity to 
the NJL model we will restrict to the structure of the phase diagram only. Following Refs. 0, 29 1 
we choose the model parameters by /jr = 93MeV, Mg = SOOMeV, rn-o- = 600MeV and vary the 
pion mass m-,^. In Figs. I5|6I we show the phase diagrams including the domain for inhomogeneous 
phases with m.,^ = 0, 69, 138MeV. The model in the applied approximation does not have a critical 
point in the chiral limit and as expected after the discussion of the model in section IVIII and the 
numerical results for the NJL model, the first order phase transition from the chirally broken to 
the (partially) restored phase is replaced by two second order phase transitions at the boundary 
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of a inhomogeneous domain. For half the physical pion mass m.,^ = 69MeV we find the same 
scenario, whereas we find a critical point in the phase diagram at the physical value of the pion 
mass = 138MeV. Although the critical point is in principle hidden by inhomogeneous phases 
as discussed at the end of section IVIH the Lifshitz point at the cusp of the two second order 
phase transition lines and the critical point appear to be so close that we cannot distinguish those 
numerically. 

IX. DISCUSSION 

In this paper we have shown how inhomogeneous phases in the NJL and the QM model with lower 
dimensional modulations can be investigated and especially how this analysis is related that in the 
(chiral) Gross-Neveu model. As suggested in Ref. [13] we have confirmed that the first order phase 
transition line in the NJL phase diagram of homogeneous phases is completely replaced by two 
second order phase transition lines that border an inhomogeneous phase. The intersection of these 
lines defines a Lifshitz point which coincides with the critical point also for finite current quark 
masses. For the QM model in the applied approximation this picture is somewhat modified as the 
critical point (in the case of homogeneous phases) is slightly hidden by an inhomogeneous phase. 
Although the qualitative behavior is very general, the specific structure of the phase diagram and in 
particular the importance of inhomogeneous phases is strongly depending on the model parameters. 
We have established that there is a close connection between the significance of inhomogeneous 
phases and the strength of the first order phase transition when limiting to homogeneous phases, in 
particular the size of its spinodal region. In addition also the relevance of pseudo-scalar condensates 
has been addressed. 

This work is a further step in the investigation of inhomogeneous phases within NJL-type models. 
In the following various related problems could be addressed: From a phenomenological point of 
view it would be interesting to study the interplay and competition with color-superconducting 
phases, being expected in the same domain of the phase diagram [l|, 0]. Also the calculation of 
higher-dimensional modulations should be in range of numerical approaches, especially since the 
magnitude of the order parameter and the size of the Brillouin zone can be large and a numerical 
determination of the density of states therefore be possible [30]. More towards an extension of 
the considered models, a GL analyses including e.g. vector-vector and 't Hooft interaction could 
be performed to check whether and how the relation of Lifshitz and chiral critical point can be 
modified. Also inhomogeneous phases in the Polyakov-NJL model can be studied easily since the 
density of states for quasi-particles is not altered by the inclusion of the Polyakov loop. As a con- 
sequence the technical issues associated with inhomogeneous phases are not furtherly complicated 
here. From a technical point of view e.g. the zero point energies in the QM model can be easily 
taken into account since the asymptotic behavior of the spectral density is known analytically. This 
would then allow for a proper renormalization of the model. In general it would of course also be 
interesting to study these phases beyond mean-field approximation, especially since fluctuations 
play a key role at a critical point. 
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